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Abstract 

This paper is devoted to the study of solitary waves and solitons whose 
existence is related to the ratio energy/charge. These solitary waves are 
called hylomorphic. This class includes the Q-balls, which are spherically 
symmetric solutions of the nonlinear Klein-Gordon equation (NKG), as 
well as solitary waves and vortices which occur, by the same mechanism, 
in the nonlinear Schroedinger equation and in gauge theories. This paper 
is devoted to the study of hylomorphic soliton. Mainly we will be inter- 
ested in the very general principles which are at the base of their existence 
such as the Variational Principle, the Invariance Principle, the Noether 
theorem, the Hamilton- Jacobi theory etc. We give a general definition 
of hylomorphic solitons and an interpretation of their nature (swarm in- 
terpretation) which is very helpful in understanding their behavior. We 
apply these ideas to the Nonlinear Schroedinger Equation (NS) and to the 
Nonlinear Klein-Gordon Equation (NKG) repectively. 
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1 Introduction 

Roughly speaking a solitary wave is a solution of a held equation whose energy 
travels as a localized packet and which preserves this localization in time. A 
soliton is a solitary wave which exhibits some strong form of stability so that it 
has a particle-like behavior. 

Today, we know (at least) three mechanism which might produce solitary 
waves and solitons: 



• Complete integrability; e.g. Kortewg-de Vries equation 

u t + u xxx + 6uu x = 

• Topological constrains: e.g. Sine-Gordon equation 

u t t - u xx + sintt = 

• Ratio energy/charge: e.g. the following nonlinear Klein-Gordon equation 



— Aip + — J— = 0; V G 
l + \ip\ 



This paper is devoted to the third type of solitons which will be called 
hylomorphic solitons. This class of solitons includes the g-balls (see [20 ) which 
are spherically symmetric solutions of as well as solitary waves and vortices 
which occur, by the same mechanism, in the nonlinear Schroedinger equation 
and in gauge theories. 

The most general equations for which it is possible have hylomorphic solitons 
need to have the following features: 

• A-l. The equations are variational namely they are the Euler- Lagrange 
equation relative to a Lagrangian density C. 
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• A-2. The equations are invariant for time translations, namely C does 
not depend explicitly on t. 

• A-3. The equations are invariant for a gauge action, namely C does 
not depend explicitly on the phase of the field ^ which is supposed to be 
complex valued (or at lest to have some complex valued component). 

By Noether theorem A-l and A-2 guarantee the conservation of energy, 
while A-l and A-3 guarantee the conservation of an other integral of motion 
which we call hylenic charge. 

The existence of hylomorphic solitons is guaranteed by the interplay between 
energy and hylenic charge. 

This paper is devoted to the study of hylomorphic soliton. Mainly we will be 
interested in the very general principles which are at the base of their existence 
such as the Variational Principle, the Invariance Principle, the Noether theorem, 
the Hamilton- Jacobi theory etc. These principles will be discussed in section [2] 
We will also give a proof of Noether theorem in the form needed for this study. 

In section [31 we give a general definition of hylomorphic solitons and we give 
an interpretation of their nature (swarm interpretation) which is very helpful 
in understanding their behavior. 

In sections H] and El we apply these ideas to the Nonlinear Schrocdingcr Equa- 
tion (NS) and to the Nonlinear Klein-Gordon Equation (NKG) repectively. The 
amount of results relative to NSE end NKG is huge (we just to cite some of the 
main papers: [35], [33], [37], QJ], [H], [IS], [36], [2D] and the books [35], [33])); 
however, here we restrict ourselves to the aspects of these equations related to 
the ideas of the previous sections; in particular, we will emphasize their com- 
mon "hylomorphic" nature which has been recently discovered. The technical 
aspects related to the study of these equations will be treated superficially, but 
we will refer the interested readers to the appropriate papers. 

An other set of equations to which this kind of ideas can been applied are 
the Klein-Gordon-Maxwell Equations (NKGM) and also more general gauge 
theories such as the Yang-Mills equations. In this paper, we will not discuss 
these equation and we refer to [5J, [7J, [16], [23], [24], [32], [9], [10] and their 
references. 

An other very interesting aspect of NS, NKG and NKGM is the fact that 
they admit hylomorphic vortices namely solitary waves having a nonvanishing 
angular momentum. The existence of these vortices is guarateed by a mechanism 
similar to that of hylomorphic solitary waves. Here we will not discuss the 
hylomorphic vortices and we refer to [13], [5], [5], [5] and [ID] . 

2 The general theory 
2.1 The variational principle 

The fundamental equation of Physics are the Euler-Lagrange equations of a 
suitable functional. This fact is quite surprising. There is no logical reason for 
this. It is just an empirical fact: all the fundamental equations which have been 
discovered until now derive from a variational principle. 

For example, the equations of motion of k particles whose positions at time 
t are given by Xj(t), Xj S R 3 , j = 1, k are obtained as the Euler Lagrange 
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equations relative to the following functional 




where rrij is the mass of the j-th particle and V is the potential energy of the 
system. 

More generally, the equations of motion of a finite dimensional system whose 
generalized coordinates are qj(t), j = 1, k are obtained as the Euler Lagrange 
equations relative to the following functional 



S = J >qk,qx, ,4k) dt 



Also the Dynamics of fields can be determined by the Variational Principle. 
From a mathematical point of view a field is a function^ 

u : R Ar+1 — > R fc , u={ux,....,u k ). 

where R JV+1 is the space-time continuum and R fc is called the internal parame- 
ters space. Of course, in physical problems, the space dimension N is 1, 2 or 3. 
The space and time coordinates will be respectively denoted by x — (x%, ..,xn) 
and t respectively. The function u(t, x) describes the internal state of the ether 
(or vacuum) at the point x and time t. 

From a mathematical point of view, assumption A-l states that the field 
equations are obtained by the variation of the action functional defined as fol- 
lows: 

S= C(t,x,u,Vu,d t u) dxdt. (2.2) 



The function C is called Lagrangian density function but in the following, as 
usual, we will call it just Lagrangian function. 

If u is a scalar function, the variation of (|2.2p gives the following equation: 



A d ( dC \ dC 



IT. = (2.3) 



^— ' dxi V du x . / du 

If u = (ui, ....,Uk), the Euler-Lagrange equations take the same form provided 
that we have use the convention that 

dC ( dC dC \ dC ( DC DC 



du Xi \du\ tX ^ ' duk, Xi J ' du \dui ' duk, 
So, if u has fc components ( k > 1) then eq. (|2.3|) is equivalent to the k equations: 



V— ( dC \ - — 

^ dxi \ du Lxi j due 



0, £=l,...,fe (2.4) 



use the convention to use Greek letters ip, ^/ etc. to denote complex valued functions 
and Latin letters u, v,... to denote real valued function. 
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2.2 The invariance principle 

A functional J is called invariant for a representation T g of a Lie group if 

J(T g u) = J{u). (2.5) 
Now, let us consider the variational equation 

F{u) = 

where F(u) = dJ{u). If J is invariant, given any solution u, we have that also 
T g u is a solution. 

We need to be careful in the interpretation of Q2.5jl . In fact, if u belongs to 
some function space $ (ft, V), (where ft C R Ar+1 and V is a finite dimensional 
vector space), it might happen that T g u $ {ft, V). For example, if 

(T h u) (x) = u(x -h), he R N+1 

we have that T^u G $ (ft — h, V) . 

Thus we are led to give the following definition: we say the Lagrangian 
C (t, x, u, Vw, dtu) is invariant with respect to the representation T g if 

C (t, x, u, Vu, d t u) = C (t', x', vf, Vu', d t u') (2.6) 

where u' (t' , x') — T g u (i, x) . 

In this case, the equation (|2.5|) need to be interpreted as follows: 



/ £(t',x',u',\/u',d t u')dxdt= / C{t,x,u,^u,d t u) dxdt (2.7) 
JT g n Jn 

where 

T g n := {{t',x') G W N+1 : (t,x) G fl} . 

We say that a Lagrangian L is invariant for the action T g if (|2.7[) holds for all 
bonded sets ft C R N . 



2.3 The Poincare invariance 



The fundamental equations of Physics are invariant for the Poincare group: it 
is the basic principle on which the special theory of relativity is based. 

The Poincare group *p is a generalization of the isometry group l£. The 
isometry group 2; in is the group of transformation which preserves the 
quadratic form 

N 
i=l 

i.e. the Euclidean norm and hence the Euclidean distance 



iV 



\ ^2\xi-yiY 
\ i=l 



namely, if g G (£, 



dE(x,y) 

d>E(gx,gy) = d E (x,y). 
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If we identify the physical space with R 3 , the isometry group is also called 
congruence group. Roughly speaking, the content of Euclidean geometry is 
the study of the properties of geometric objects which are preserved by the 
congruence group. 

The Poincare group *P, by definition, is the transformation group in R^ 1 
which preserves the quadratic form 

2 



i=i 

which is induced by the Minkowski bilinear form 

N 

(x, y) M = ~ x oyo + XiVl 



The Minkowski vectors v = (v , ..,%) = («o, v ) are classified according to 
their causal nature as follows: 

• a vector is called space-like if (v, v) M > 0, 

• a vector is called spacedike if light dike if (v, v) M = 0, 

• a vector is called space-like if time-like if (v, v) M < 0. 

The casual nature is not changed by a Poincare transformation, and hence 
it is not a transitive group (as the isometry group): space and time-are mixed, 
but not so much. 

In the real world we have N = 3 and the Poincare group is a 10 parameter 
Lie group generated by the following one-parameter transformations: 

• Space translations in the directions x,y,z : 

x' = x + x x' = x x' = x 

y' = y y' = y + yo . y' = y 

z' = z ' z' = z ' z' = z + zq 
t' = t t' = t t' = t. 

This invariance guarantees that space is homogeneous, namely that the 
laws of physics are independent of space: if an experiment is performed 
here or there, it gives the same results. 

• Space rotations: 

x' = x x' = x cos 6 2 — z sin 8 2 x' = x cos 9 3 — y sin 9 3 

y' = y cos 9i — z sin Q\ y' = y y' = x sin 8 3 + y cos 63 

z' = y sin Q\ + z cos 9\ ' z' = x sin 9 2 + z cos 9 2 ' z' = z 

t' = t t' = t t' = t. 

This invariance guarantees that space is isotropic, namely that the laws 
of physics are independent of orientation. 
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• Time translations: 

x' = x 

?/ = y 
z' = z 

t' = t + t . 

This invariance guarantees that time is isotropic; namely that the laws of 
physics are independent of time: if an experiment is performed earlier or 
later, it gives the same results. 

• Lorentz boosts: 

x' = 7 [x — v±t) x' = x x' = x 

y' = y . y' = i(y- v 2 t) . y' = y 

z' = z ' z' = z ' z' = 7 (z — v^t) 

t' = j(t-vix) t' = j(t-v 2 y) t' = j(t-v 3 z) 

where 



(2. 



-"Trh- (2 - 9) 

with v = Vi, i — 1,2,3. This invariance is an empirical fact and, as it will be 
shown in section 15. 5( it implies the remarkable facts of the theory of relativity 
such as the space contraction, the time dilation and the equality between mass 
and energy.. 

The Lorentz group is the 6 parameters Lie group generated by the space 
rotations and the lorentz boosts (plus the time inversion, t — > — f, and the 
parity inversion (x,y,z) — > (—x,—y,—z)). Clearly it is a linear subgroup of 
GL(6). 

The Poincare group is the 10 parameters Lie group generated by the Lorentz 
group and the space-time translations. Then it is a subgroup of the affine group 
in M 6 . 

The Poicare group acts on a scalar field by the following representation: 

(T g iP)(t,x)=ijj(t',x'), (t',x')=g(t,x) 

The simplest equation invariant for this representation of the Poincare group 
is the D'Alembert equation: 

= (2.10) 

where 

<9 2 V> a i j a , d 2 i> d 2 ^ <9V 

U ^ = -t^t ~ A ^ and A ^ = ttt + ttt + in:- 

ot A ox z ay 1 oz z 

The D'Alembert equation is the simplest variational field equation invariant 
for the Poincare group. 

In fact it is obtained from the variation of the action 

S v = -\J (#.#>m dxdt = ^J [\d t ij\ 2 - |V^| 2 ] dxdt. (2.11) 

In this case, the Lagrangian function is given by 

A) = -~ (#, #> M = \ IftVI 2 - \ |VV| 2 (2.12) 

It is easy to check that if ip is a solution of this equation, then also T g ip is a 
solutions of the equation for every g £ 
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2.4 The Galileo invariance 

The Galileo group (3 as the Poincare group is a trasformation group on the space 
time R 4 . The Galileo group, by definition, if the set of trasformations which 
preserves the time intevals and the Euclidean distance between simultaneous 
points. More pricisely, an affine trasformation g s © if, given two points (ti, xi) 
and (t2, X2) , we have that 

t'l — t' 2 = ti — t'2 

and 

ti = t% (ii = t 2 and d,E {x'i — x' 2 ) = A,e (x\ — X2)) 

where (t^x^) = g(ti,X{) and dE{x,y) is the Euclidean distance. 

Thus the Galileo group is a 10 parameters Lie group generated by the space- 
time translation the space rotations but the Lorentz boosts (|2.8p are replaced 
be the Galilean tranformations namely by the tranformations 



x' — x — v\t x' — x x' — x 

y' = y . y' = y - v 2 t . y' = y 

z' = z ' z' = z ' z' = z — v$t 

t' = t t' = t t' = t 



(2.13) 



The equations of classical mechanics are invariant for the Galileo group. We 
are interested in field equations which are invariant for a representation of the 
Galileo group. 

Given the Galileo tranformation g v : R 4 — > R 4 defined by 

g v (t,x) = (t,x- vt). 
we consider the representation T gv : L 2 (W 4 , C) — > L 2 (R. 4 , C) defined by 

(T 9 » (t, x)=1>(t,x- vt) e i( y x -^ H \ (2.14) 
The Shroedinger equation for a free particle (of mass 1) 

is the simplest second order equation invariant for a representation of the Galileo 
group on the space of complex valued vector fields L 2 (M 4 ,C). 

2.5 The Gauge invariance 

Take a function 

tp : R 4 -> V 

and assume that on V acts the representation T g of some group (G,o) . This 
action induces two possible action on ip : 

• a global action: tp (x) 1— > T g ip (x) where g E G 

• a local action: (x) 1— > T g ^ip (x) where g(a;) is a smooth function with 
values in G. 
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In the second case, we have a representation of the infinite dimensional group 

= C°° (R 4 ,G) 

equipped with the group operation 

(goh) (x) = g(x) o h(x) 

If a Lagrangian C satisfies the following condition, 

C(t, x, rj), VV, W) = C (t, x, T ff </>, V (T g i/>) , d t (T g ip)) , geG 

we say that it is invariant for a local action of the group G, or for a trivial gauge 
action of the group G; if C satisfies the following condition, 

C(t, x, ip, Vi/>, d t ip) = C (i, x, T g ( x )ip, V (T g{x) ip) , d t (T g(x) ip)) , g(x) e<S 

we say that it is invariant for a local action of the group G, or for a gauge action 
of the group 0. 

Let us consider two simple examples: the functional 



J C (Vu) dx , u e 



is invariant for a global action of the group (M, +) . In fact, if we set T r u = u + r, 
r £ R, we have that 

£{Vu) =£(V (T r u)). 
Next, consider the functional 

C (da) dx 

where a is a 1-form and d is the exterior derivative of a. In this case, C (da)is 
not only invariant for a trivial action of (R, +) , but also for the local action 

T g{x) a = a + dg{x), g (x) e © := C°° (M 4 , R) 

in fact 

C(d(a + dg{x))) =C(da) . 

The simplest gauge invariance can be obtained taking a complex valued 
scalar field 

ip : R 4 -> C, 

and to consider the group S 1 = {e i6 : 9 £ R} and the following representation 

i/>^e l V (2-15) 

The Scroedinger equation and the Klein-Gordon equation are invariant for the 
global action (|2.15|) . The Klein-Gordon-Maxwell equations are invariant for a 
local action (|2.15[) . For a discussion of these aspects of KGM, we refer to [?J]. 
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2.6 Noether's theorem 



In this section we will give a proof of Noether theorem stated in a suitable form 
for the applications considered in this paper. 
First of all, we need the following lemma 

Lemma 1. Let p : R N+1 — > K and J :R N+1 -> M w 6e two smoot/i functions 
defined on the "space-time". Assume that they satisfy the continuity equation 



and that for all t 



Then for all t 



Proof. Let 



dp, 



p{;t), wM) andJ(-,t) are in L 1 (M 3 ) 



d f 

— / p(x,t)dx = 
at j r n 



(2.16) 

(2.17) 
(2.18) 



B R = {x e M 3 : \x\ < R} , R>0 
then, integrating on Br, we get 
dp 



dt 



dx = — / V • Jdx = — 



[ (J • n)da 



(2.19) 



where n denotes the outward normal to the boundary OBr of Br. Then 

| J • n| da (2.20) 



b r dt 



< 



f)B r , 



Since §f(.,*) is in L 1 (E 3 ), there exists lim J Br ^dx 
prove that this limit is 0. Arguing by contradiction we assume that 

dp 



lim 

Ft—. >oo 



at 



dx 



= a > 



and we have to 



(2.21) 



then, by (f2T20|) and (|2~2T]l . the map ip defined by 



<p(R)-- 

is not integrable in (0, +oo) and 

| J ■ n| dx = 

which contradicts assumption (|2.17|) . 
□ 



|J • n| da 



8B n 



+ 00 



ip(R)dR = +oo 



Suppose that a Lagrangian is invariant for the action T g of some Lie group G. 
We denote by T s (a) € K) the action of a one-parameter subgroup {g (A)} AgR . 
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Notice that this subgroup is isomorphic either to S 1 or R. We use the following 
notation: 

u x = T g(x) u (2.22) 
and, if the group acts also on the independent variables, we set 

t x = T gW t (2.23) 
xx = T g{x) x (2.24) 

For example, consider the hrst of the Lorentz transformation (|2.8[k in this 
case the parameter A is the first component of the velocity v, and we have 

t — vx 

tv = 



X\ — vt 



VT^2 

u v = u (t v ,x l7V ,x 2 ,x 3 ) 

Lemma 2. If £ is invariant with respect to a one parameter group g (A) , then, 
using the notation S2.22\) . V2.23\) . \2.2J$ , we have that 



_d_ 



£(tx,xx,ux,Vux,d t ux)<fi(tx,xx) dxdt = 



where tt = [t , t x ] x R N and tp € V (0)0. 

Proof. We approximate tp with a function tp £ defined by 

3 

where xsij is the characteristic function of 

% := {(t,x) G R N+1 : je < tp(t,x) < (j + 1) e} 
If A is small, so that the support of tp (tx, xx) is contained in f2, we have that 

/ C(tx,xx,ux,^ux,d t ux)tps (tx,x\) dxdt = y^je / £ (t\,x\,u\,\/ux,d t u\) dxdt 
Jn j JT^sij 

and hence, by (|2.7[) . 

je / C(tx,xx,ux,Vux,dtUx) dxdt = /Jje / £(t,x,u,\7u,dtu) dxdt 

£ (t, x, u, Vit, dtu) tp e (t, x) dxdt 



and so, 



£ (tx, xx,ux,Vux, dtux) (tx, xx) dxdt = / £ (t, x,u,\7u,dtu) tp (t, x) dxdt 
a Jn 



2 T> (Q) denotes the space of infinitely diffcrcntiablc functions with compact support in Q. 
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Taking the limit for e — ► 0, we get that 
C (t\, x x ,u Xl Vu A , d t u x ) <p (t\, x\) dxdt 



C (t, x, u, Vu, dtu) ip (t, x) dxdt 



and we get the conclusion. 

□ 



The above lemma is useful since the introduction of the compact support 
function ip allows us to work on a fixed domain Q and hence we do not have to 
consider the variation of the domain T\£l. 

Theorem 3. Let C be invariant with respect to a one parameter group g (A) , 
and let u = u\ be a smooth solution ot the Euler- Lagrange equation h2.S\) . Using 
the notation QUfflj) , (E2$ , we set 



- f — dux _ r — 
P ~ \duxt OX dX 



(2.25) 



A=0 



N 



dC du x _ dx\ 
du x , x i dX dX 



A=0 



Then we get the continuity equation 

dp 
dt 



V • J =0. 



(2.26) 



(2.27) 



Proof. By the invariance of the Lagrangian and lemma [21 we have: 
J C<p dxdt = \/ip G T> (ft) 



d 

dX 



and hence 



' X ip + C^- I dxdt = V^eD(!l) 



dX 



dX 



(2.28) 



Let us compute each derivative; in this computation, we write u, x,t instead of 
u\,x\,t\ to make it readable and we use the notation x° = t: 



N 



dC dC d 2 u 

dX = ^ 



i=0 
N 

= E 

i=0 
N 

=E 

Then, by equation 



dC du 



du T i dXdx 1 du dX 



' d 


f dC 


du\ 


d 


( dC x 


\ du 


dx l 


\du x i 


dXJ 


dx l 


\du x i / 


/ dX 



dC du 
du dX 



d ( dC du 



dx % V du~i dX 



' N 

E 

,i=0 



d ( dC 



dx % \du x i 



dC 
du 



du 
dX 



dC 
dX 



N 

E 

i=0 



d ( dC du 



dx % V du rr i dX 



(2.29) 
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Also 



JV 



dip ^ d<P 9x l 



dx 1 dX 



Then, by (j2.28[> . (|2.29p and the above equality, we have that 



' JV 

E 



dC d 



d 



dx 1 \ du, r i dX 



N 



£ E 



dip dx % 
dx l dX 



dxdt = 



(2.30) 



Li=0 v ' i=0 

Moreover, since ip has compact support, by the divergence theorem we have that 

JV 



d 



so 



JV 

E 

i=0 



d_ 

dx 1 



C 



i=0 



dx 1 
~dX 



,dx 



ip dxdt 



JV 



c^^^d. 



By this equality and (|2.30p . we get that 

JV 







y- 

i=0 
N a 

y— 



dC du 

dx 1 V du x i dX 

dC du 



i=0 



N 



dx 1 dX 



i=0 



dx 1 V du. r i dX 



<p dxdt 

_ r— 
dX 



i=0 



dxdt 



dx 1 



E-—^ ( C—rr- ) tpdxdt 
dx l \ r)\ I r 



OX 



By the arbitrariness of ip we get 



E 



d f 8C du ^dx % 
dx 1 \ du^i dX dX 



= 



or, 



d ( dC du 



dt 



dt\du~ t dX~ 



dx 



N 



d ( dC du 



dx 1 \ du T i dX 



,dx % 
'~dX 



Then the functions (|2.25p and ()2.26[) satisfy the continuity equation (12.27p . 

□ 

Then by lemma [T] and Th. [31 we have the following result: 

Theorem 4. (Noether's theoerem) Let C be invariant with respect to a one 
parameter group g (A) , and let u be a smooth solution of the Euler- Lagrange 
eguation 112.3$ . Suppose that u decays sufficiently fast so that \2.11ty holds. 
Then, using the notation 12.22$ . 12.29$ . {2.24% we have that 



l{u) = J 
is an integral of motion. 



dC du x _ dt x 
du\ t dX dX 



dx 



A=0 
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2.7 Conservation laws 



Now, using Noether theorem [H we can compute the main integral of motion. 
They are due to the homogeneity of time and the homogeneity and isotropy of 
space which provide the invariance with respect to the time translations, space 
translations and space rotations. We consider the case in which C depends 
on a complex valued scalar function ip. The computation can be done setting 
ip = Ui + iv-2, and considering £ as function of u = (u\, U2) ■ 

• Energy. Energy, by definition, is the quantity which is preserved by the 
time invariance of the Lagrangian; it has the following form 

£ = RC J {djk)^- C ) dX (2 ' 31) 

• Momentum. Momentum, by definition, is the quantity which is pre- 
served by the space invariance of the Lagrangian; the invariance for trans- 
lations in the Xi direction gives the following invariant: 

dC 



Pi = - Re / • d Xi ipdx 

J d(d t ip) 

The numbers Pi are the components of the vector 

P = -Re/— .V^* (2.32) 

Angular momentum. The angular momentum, by definition, is the 
quantity which is preserved by virtue of the invariance under space rota- 
tions of the Lagrangian with respect to the origin 



M=Re j im)- {xx ^ )dx - 



Proof. First, we compute M3. Setting x = (x,y,z), we have that 

x\ = x cos A — y sin A 
yx = x sin A + y cos A 

z\ = z 
t\ = t 

then, setting ?/; — ui + IU2, 

dC du\ t dC du2 1 „dt\ _ / dC dip „dt 

- C— ] = Re — ; — — — L 



dC ( dip dip 



dux t t d\ du2,t d\ d\J X=Q ' \dip t d\ d\ / A=0 

= Re 

Analogously, we have 



d (d t ip) V dx ^ dy X 



\l Re / dC ( d ^ 

d {dtip) \ dz dy 



dC (dip dip 

M2 = Rel d(M)'\^ z ~^ x 
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Then we get the conclusion. 

□ 



Hylenic charge. The hylenic charge, by definition, is the quantity which 
is preserved by by the trivial gauge action (I2.15p . The charge has the following 
expression 

n = lm Jdm^ dx (2 - 33) 

Proof. We have that ip\ = tpe tX ; then 



= Im {am* 



□ 



2.8 The Hamilton-Jacobi theory 

In order to understand the motion of hylomorphic solitary waves, it is necessary 
to know the basic notions of the Hamilton-Jacobi formulation of the laws of 
Mechanics. Then, in this section, we will briefly recall these notions. 

The Lagrangian formulation of the laws of the Mechanics assumes a function 

C = C(t,q,q) 

of the generalized coordinates of the system q = (<7i, of their derivatives 
q = (qi,...,qk) and of time. The trajectories q(t) such that q (to) — xq and 
q (ti) — x\ are the critical points of the action functional 

S(q)= ^ L(t,q,q) dt (2.34) 

defined on the space 

Cx„,x! [*Q)*i] = {<? e C 1 [to,*i] : q(to) = xn and q(tx) = x x ) 
Thus a trajectory q(t) satisfies the "Euler-Lagrange" equations: 

1| -I — 1 * <"•> 

However, this is not the only formulation of the lows of Mechanics. An other 
very important formulation can be obtained as a first order system provided that 



dq 



2 >0. (2.36) 



In this case, we set 

dC 

Pj = gT- (t, q, q) (2.37) 
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By (|2.36p . we have that the function 



dqj 

is smoothly invertible and hence there exits a smooth function F such that 
(|2.37|) can be rewritten as follows: 

q = F(t,q,p) (2.38) 

Now, we can define the Hamiltonian function as follows: 

H{t,q,p) = [{p,q) - C(q,q,t)]. =F{p) 

where 

fc 

(P,q) =^2PjQj 
j'=i 

denotes the paring between the tangent space (to the space of the q's) and the 
relative cotangent space. 

Then, the action (|2.34p can be rewritten as follows 



S(p,q)= I' [(p,q) -H{t,q,p)] dt 



and the relative " Euler-Lagrange" equations take the form 

. m 

1 = ^(t,q,p) 

on 

p = -—(t,q,p) 

This is the Hamiltonian formulation of the laws of Mechanics and the above 
equations are called Hamilton equations. An other equivalent formulation of 
the laws of Dynamics is given by the Hamilton- Jacobi theory which uses notions 
both of the Lagrangian formulation and the Hamiltonian one. It reduces the 
laws of Mechanics to a partial differential equation and to a first order ordinary 
differential equation. The Hamilton-Jacobi theory has been very useful to relate 
the laws of Optics to Dynamics. For us, it essential if we want to understand 
the motion of solitons regarded as material particles. 

The starting point is the definition of a function S = S(t,x) called action. 
We fix once for ever a point (to, xq) and a point (t, x) which will be considered 
variable. Moreover, we set 

S(t,x) = f £{t,q x ,q x ) dt (2.39) 

where q x is a critical point of (|2.34p on the space C XQ x [to, t] ; in general, this point 
is not unique; however, if (|2.36|) holds, it is possible to prove the uniqueness of 
the minimum provided that (t, x) is sufficiently close to (to, Xo) {t ^ to). Hence, 
there exists an open set fi in which the function (12.39)1 is well defined. The 
function (|2.39p is called action as the functional (|2.34p . However, even if they 
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have formal similar definitions, they are quite different objects: S is a functions 
of k + 1 variables defined in an open set f2 C R^ 1 while S is a functional 
defined in the function space C]. Q X [to, t]. 

The Hamilton- Jacobi theory states that the function S, in O, satisfies the 
following partial differential equation: 

d t S + H (t, x, VS) = (2.40) 

Moreover, this result can be inverted in the sense stated by the following 
theorem: 

Theorem 5. Let S be a solution of eq. \2.4-Uty in and let q = (q±, qu) be a 
solution of the following Cauchy problem for (t, x) G Q: 

dC 

— (t,q j ,q j ) = VS(t,q j ), j = l,...,k (2.41) 
°1j 

qit) = x 

with (t, x) € O. Then, q satisfies eq. \2.35\) with initial conditions 

q(t) = x 

q(t) = F(i,x,VS(i,x)) 

where F is given by 12.38]) . 

The proof of this theorem can be found in any book of Classical Mechanics; 
for example in the beautiful book of Landau and Lifchitz 30J . 

Notice that the above Cauchy problem is well posed, at least for small times, 
since, by ((235)) . eq. (|2~4T|) gets the form 

q = F(t,q,VS(t,q)) 

where F is the smooth function given by (|2.38p . 

Thus, we can say that the equation of motions (|2.35p are equivalent to the 
set of equations, 

d t S + H(t,x, VS) = (2.42) 

q = F(t,q,VS(t,q)) (2.43) 

If C does not depend on t, then Ti. is a constant of motion (namely it is the 
energy of the system). In this case, by eq. (|2.40|) . d t S — —h, namely it does 
not depend on time and it represents the energy of the system with the sign 
changed. In this case, eq. (|2.40|) takes the form 

H (x, VS) = h. 

Let see some examples: 
Newtonian dynamics: 

£(*>(?)?) = ~ v (i) 

Then 

p = mq 
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™™ n= ^ p2+v{q) 

and equations (12.42|2.43p take the form 

d t S+^-\VS\ 2 + V(x) =0 (2.44) 
2m 

q=-VS(t,q) (2.45) 
m 

Relativistic dynamics: The Lagrangian of a relativistic particle is given 

by: 

£ (t, q, q) = -m \/l - Q 2 

where mo is a parameter. We refer to Landau-Lifchitz [31] for a very elegant 
deduction of this Lagrangian from the Mikowsky geometry of space-time. 
Then 

d£ m . . 
P = -7TT = 7= q = jm q 2.46 
dq y/l - q 2 

with 

1 

7 



1 

v A y 

and equation (|2.38|) becomes 



« = / fe 2 ( 2 - 4? ) 

Since the mass of a particle is defined by the equation m — p/q, we will get 
that the mass changes with velocity 

m = 7mo 

and the interpretation of mo as rest mass. The Hamiltonian is: 

v ■ i R ^ TO ° -2 , /i ^ ™o 

H=pq + m vi - g — , „ g + moV 1 - <r = > = 7 TO o 

v 1 -q v 1 - ? 

Since the Lagrangian in independent of time, the Hamiltonian represents the 
energy and this gives the Einstein equation: 

£ = H = m = 7m (2.48) 

Now let express H as function of p. Using cq. (12.47j) we get 



1 - 3^—7 



H (p, 9) = . m ° = = t/mg + P 2 (2-49) 



and the equations (|2.42I2.43| take the form 



flfcS + ^mg + |VS| =0 (2.50) 
• ^ (2.5D 
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3 Hylomorphic solitary waves and solitons 



3.1 An abstract definition of solitary waves and solitons 

Solitary waves and solitons are particular states of a dynamical system described 
by one or more partial differential equations. Thus, we assume that the states 
of this system are described by one or more fields which mathematically are 
represented by functions 

* : R N -» V (3.1) 

where V is a vector space with norm \ ■ \ v which is called the internal parameters 
space. We assume the system to be deterministic; this means that it can be 
described as a dynamical system (X, U) where X is the set of the states and 
U : R x X — > X is the time evolution map. If ^o(x) € X, the evolution of the 
system will be described by the function 

*(t,x) = U t *o(x) (3.2) 

Now we can give a formal definition of solitary wave: 

Definition 6. A state "Jo G X, is called solitary wave if its evolution has the 
following form: 

where ht and gt are transformations of V and M. N respectively. 

For example, consider a solution of a field equation which has the following 
form 

(t, x) = e-^-^^oix - vt); * G L 2 (WL N , C); 

then the conditions of the above definition are satisfied with h t ip — e _l ( k ' x_u '*V 
(tp € C) and g t x = x — vt. 

The solitons are solitary waves characterized by some form of stability. To 
define them at this level of abstractness, we need to recall some well known 
notion in the theory of dynamical systems. 

Definition 7. Let X be a metric space and let (X, U) be a dynamical system. 
An invariant set T C X is called stable, if Ve > 0, 3S > 0, V* e X, 

d(*,r)<s, 

implies that 

Vt e M, d(y t V,T) < e 



Definition 8. A state \Po is called orbitally stable if there exists a finite di- 
mensional invariant stable manifold V such that G T- 

The above definition needs some explanation. Since T is invariant, Ut^o G 
r for every time. Thus, since T is finite dimensional, the evolution of is 
described by a finite number of parameters. Thus the dynamical system (r, U) 
behaves as a point in a finite dimensional phase space. By the stability of 
T, a small perturbation of remains close to V. However, in this case, its 
evolution depends on an infinite number of parameters. Then, as time goes 
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on, the evolution of the perturbed system might become very different from 
Ut^a- Thus, this system appears as a finite dimensional system with a small 
perturbation. We refer to section (|4. 7[) where this fact will be seen in details in 
a concrete case. 

Definition 9. A state \l/o € A, / 0, is called soliton if it is a orbitally 
stable solitary wave. 

In general, dim (r) > N and hence, the "state" of a soliton is described by 
N parameters which define its position and other parameters which define its 
" internal state" . 

3.2 Definition of hylomorphic solitons 

Now let us assume that our system satisfies assumptions A-l, A-2 and A-3 
of pagf2] By Noether theorem, assumptions A-l and A-2 of guarantee the 
conservation of the energy £ ('J) (see (|2.31jl ) while A-l and A-3 guarantee the 
conservation ot the hylenic charge Tt ('I') (see (|2.33[) ). 

Definition 10. A stationary wave 'J'o is called hylomorphic wave if 

£ (1r ) = c„ := min £ (*) for some a e K (3.3) 

where 

Moreover <3>o is called hylomorphic soliton if it satisfies definition [5] namely if 

T a = {V eX:£ (*)=c ff } (3.4) 
is a finite dimensional stable manifold. 

Observe that the energy and the charge are constant of the motion, then T a 
is an invariant set. Then the above definition is consistent with definition [9l 
Let <J/o be a hylomorphic stationary wave as in Def . [TU] and set 

V(*) = (£ (*)-c a ) 2 + (H(*)-*) 2 ; (3.5) 

then we have that 

Proposition 11. Let an d To- be as in Def. \1U\ and let be any sequence 
in the phase space X which is supposed a metric space with metric d. Then if 

v(y n )^o^d(v n ,r a )^o (3.6) 

^0 *s a hylomorphic soliton. 

Proof. By (|3.6p . it is immediate to check that V is a Liapunov function 
relative to the flow Ut defined by 13.21 Then, by the well known Liapunov 
theorem, it follows that T a is stable and by the definitions[S]and[n]the conclusion 
follows. 

□ 
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If is a stationary wave of a dynamical system relative to a Lagrangian 
which is invariant for the Lorentz or the Galileo group, then it is possible to 
obtain a travelling wave just making a Lorentz boost (see (|2.8I) ) or a Galilean 
transformation (see l|2.13|> ). More precisely, let T v be the representation of a 
Lorentz boost (or a Galilean transformation) relative to our system and let 

¥(t,aO = Di¥o(aO 

be the evolution of our stationary wave Wo (x) ; then 

y'(t',x') = T v *(t,x) 

is a solution of our equation which moves in time; then \I/ V (x) := *&'(t' ,x')\t=o 
is a travelling wave at the time t = 0. At pag. and pag. [33 we will see as 
this principle works in some particular cases. 

Obviously if ^q(x) is a hylomorphic standing wave \& v (x) will be called 
hylomorphic travelling wave; moreover, if ^o(x) is orbitally stable, also \I/ V (x) 
is orbitally stable and hence is a hylomorphic soliton. 

3.3 Structure of hylomorphic solitons 

Now, as it happens in Th. HJ we assume £ and H. to be local quantities, namely, 
given f 6l, there exist the density functions ps,^ (x) and pn,y (x) € L 1 (R JV ) 
such that 

S (*) = y p £> * (x) dx 

h (*) = y p-h,* (a) ^ 

Energy and hylenic density allow to define the density of binding energy as 
follows: 

/3(i,x) = #p(t,x) = [Eo ■ \ph,v (t,x)\ - ps ,* (i,x)] + (3.7) 

The support of the binding energy density is called bound matter region; 
more precisely we have the following definition 

Definition 12. Given any configuration W, we define the bound matter re- 
gion as follows 

E (*) = {x:/?*(t,x)^0}. 
If is a soliton, the set E (^o) is called support of the soliton at time t. 

Thus a hylomorphic soliton \&o consists of bound matter localized in a precise 
region of the space, namely E (Wo). This fact gives the name to this type of 
soliton from the Greek words "hyle"=" matter" and " morphe" =" form" . 



We now set 

E = lim inf £ ^\, (3.8) 

where 

X e = {W S X : Vx, |W(x)| y < e} . (3.9) 
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Now suppose that there is a state ^ which satisfies the inequality 



< En 



which will be called hylomorphy condition. The quantity 



A(¥) = 



(3.10) 



(3.11) 



is also an invariant of motion and will be called hylomorphy ratio. 

Notice that, by def. (j3.8p and p. lip , for any hylomorphic solitary wave 
we have that 

A(* ) E " 



< 

a a 



However, if (|3.10p holds, we have that 



A(*o)<^ 
a 



Actually the hylomorphy condition (|3.10p seems to be a necessary condition in 
order to have hylomorphic solitons. 

The hylomorphy condition (|3.10p guarantees the presence of bound matter 
even if no soliton is present: 

Proposition 13. If A (*(0, •)) < E Q) then for allteR 

£(*(*,.)) 



Proof. To fix the ideas, assume 7i (^f) > 0. 

P(t,x)= / [\E pH,u t * - pe,u t v (x)Y 



□ 



> / E p n ^ Ut ^ (x) - P£,u t v (x) 

= E H (E/jtf) - £ (£/ t *) 
= E H (*) - £ (*) 

e(9) 



E 



= W (*) [E - A (*)] >0 



If ^>(x) is a finite energy field usually it disperses as time goes on, namely 
\im\\U t ^(x)\\ Loo{RNy) =0 

However, if A (^) < Eq, this is not the case: 
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Proposition 14. If A($?) < E , then 

Proof: To fix the idea, set A ("5?) — Eq — a, a > 0. We argue indirectly and 
assume that, for every e > 0, there exists i such that 

II^IIl^r" v) < e 

namely, Ujty E X £ where X £ is defined by (|3.9[) . Then, by (|3.8[) . if e is sufficiently 
small 

A ( W >E - a - 

Since A (f/^) = A (if?) we get a contradiction. 
□ 

Thus if A (\&) < Eo, by the above propositions the field W and the bond 
matter field f3^ will not disperse but will form bumps of matter which eventually 
might lead to the formation of one or more hylomorphic soliton. This remark 
shows the importance of the hylomorphy condition (| 3 . 1 [) . 



3.4 The swarm interpretation of hylomprphic solitons 

Clearly the physical interpretation of hylomorphic solitons depends on the model 
which we are considering. However we can always assume a conventional inter- 
pretations which we will call swarm interpretation since the soliton is regarded 
as a swarm of particles bound together. In each particular physical situation 
this interpretation, might have or might not have any physical meaning; in any 
case it represents a pictorial way of thinking of the mathematical phenomena 
which occur. This interpretation is consistent with the names and the definitions 
given in the previous section. 

We assume that "S? is a field which describes a fluid consisting of particles; 
the particles density is given by the function pu(t,x) — pn,4>{t,x) which, of 
course satisfies a continuity equation 

dtPH + V • 3 n = (3.12) 

where 3u is the flow of particles. Hence TL is the total number of particles. 
Notice that TL does not need to be an integer number and you may assume that 
fractional particle can exist or that TL is a sort of limit valid for a very large 
number of particle as it happens in fluid-dynamics. Also, in some equations as 
for example in NKG, TL can be negative; in this case we assume the existence 
of antiparticle. 

Thus, the hylomorphy ratio 

A(,) = M 

1 ; \n(n 
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represents the average energy of each particle (or antiparticle) . The number Eq 
defined by (|3.8[) is interpreted as the rest energy of each particle when they do 
not interact with each other. If A > Eq, then the average energy of each 
particle is bigger than the rest energy; if A ('J) < Eq, the opposite occurs and 
this fact means that particles act with each other with an attractive force. 

If the particles would be at rest and they would not act on each other, their 
energy density would be 

E ■ \p n (t,x)\ ; 
but their energy density is ps(t, x); if 

ps(t,x) < E Q ■ \pn(t, x) \ ; 

then, in the point x at time t, the particles attract each other with a force 
which is stronger that the repulsive forces; this explains the name density of 
bond energy given to f3(t, x) in (|3.7p . 

Thus a soliton relative to the state 'J can be considered as a "rigid" object 
occupying the region of space E (cf. Def. [TSJ) ; it consists of particles which 
stick with each other; the energy to destroy the soliton is given by 

/ (3^(t,x)dx = / (E \p n (t,x)\ -ps(t,x))dx 

However out of E (v?) the energy density is bigger than E ■ \pn{t, x) \ ; thus 
the total energy necessary to reduce the soliton to isolated particles is give by 

E -\H (tf). 

As it is shown in Prop. [T31 there are states 'J such that E ('J) ^ but they 
are not necessarily solitons. In these states, E (^P) is a region where the particles 
stick with each other but they do not have reached a stable configurations; the 
shape of E ($) might changes with time. In may concrete situations, such states 
may evolve toward one or more solitons. In these cases we say that the solitons 
are asymptotically stable. The study of asymptotical stability is a problem quite 
involved. We refer to , [TS], [H] and [52] and their references. 



4 The nonlinear Schroedinger equation 

4.1 General features of NS 

The Schroedinger equation for a particle which moves in a potential V(x) is 
given by 

i^ = ~^ + V{x)i} 
We are interested to the nonlinear Schroedinger equation: 

i^ = -\^ + \w'{iP)+V{x)i } (4.1) 

where 

dW dW 

B "«« = Sr + C ( "» 
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namely 

w'&) = F'm^. 

for some smooth function F : [0, oo) — > K. 
We always assume that 

W(0) = W'{0) = 
If V(x) = 0, then we get the equation 

i^ = -\^ + \w'^); (NS) 

this equation can be considered as the simplest equation which is variational 
and invariant for a representation of the Galileo group. 
First of all let us check that it is variational: 

Proposition 15. Equation J^. j[ ) is the Euler- Lagrange equation relative to the 
Lagrangian density 

L = Re (idtH^) ~ \ |V^| 2 -W(ip)- V(x) \tp\ 2 (4.3) 

Proof Set 



<Si(V>) = j Re {idtipW} dxdt; 5 2 (V>) = - 
and set ip = u\ + iu 2 . We have 

S\(ip) = I Re (id t ipip) dxdt 



\\ViP\ 2 + W(^) + V(x)\iP\ 2 



dxdt 



Re [(idtUi — 9(^2) (ui — W2)] g?xg?£ 
{dtU\U2 — dtUiUx) dxdt 

Then, if ip = V\ + iv 2 

d<Si(V>) [<p] = j {d t u\V2 + d t viu 2 - d t u 2 vi - d t v 2 ui) dxdt 
(2dtU\v 2 — 2dtu 2 v\) dxdt 
Re [2i (d t ui + d t u 2 ) (v\ - iv 2 )} 
Re (2id t ipTp) 

Then 

dS 2 (ip)[tp] = - [ [Re(W,V<^) +Re(W' (ip)Tp) + 2Re(V(x)ipTp)]dxdt 
[Re (-AipTp) + Re (W (ip) Tp) + 2 Re (V(i)^)] dad* 
Rc [(-Aip + W (ip) + 2V(x)ip) Tp] dxdt. 
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dS(ip) [f \ = / Re \{2id t ^ + A^ + W' (ip) + 2V(x)ip) Xp] dxdt 
So, the critial points of S, satisy the equation 

2id t tp + Atp - W (ip) - 2V(x)ip = 



which is equivalent to (|4.1|) . 

□ 

Sometimes it is useful to write ip in polar form 

4>(t,x) = u{t, x y s{t > x) . 



(4.4) 



where u(t,x) £ R + and S(t,x) £ R/ (27rZ). Thus the state of the system ip is 
uniquely defined by the couple of variables (u,S). Using these variables, the 
action S = J Cdxdt takes the form 



S(u,S) = -J 



~ |Vu| 2 + W(u) + (d t S + ~ \WS\ 2 + V(x) ) >r 



dx 



and equation (|4.ip becomes: 



^Au+^W , (u) + (d t S+^\X7S\ 2 + V(x) )(/ = n 
ft (u 2 ) + V ■ (u 2 WS) = 



(4.5) 
(4.6) 



4.2 First integrals of NS and the hylenic ratio 

The results of section l2~7'l and easy computations show that the integral of motion 
of l4.1l are given by the following expression: 



Energy. We have 
Sty) 



^\ViP\ 2 + W(iP)+V(x)\iP\ 2 



dx 



(4.7) 



Using (|4.4|) we get: 

E(ip) = J Q |Vu| 2 + W(u)j dx + J (^\VS\ 2 + V(x) ) ird.r (-1. 

• Momentum. The momentum is constant in time if the Lagrangian is 
space-traslation invariant; this happens when V is a constant. In partic- 
ular, in eq. INS) we have that 



P = Im / Vipip dx 



(4.9) 



Using (|4.4p we get: 



P = / u 2 V5 dx 



(4.10) 
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• Angular momentum. If we assume that V is a constant, the angular 
momentum is an integral of motion and, in INSi we have 



M = Im / xx Vipip dx 



(4.11) 



Using (|4.4p we get: 



M = / xxVStt 2 dx (4.12) 
• Hylenic charge. Here the hylenic charge has the following expression 



H{ip) = / \ip\ dx = / u 2 dx 



(4.13) 



• Barycenter (or Hylecenter) velocity. The quantity preserved by the 
Galileo transformation (|2.14p is the following 



K 



i dx — tP 



(4.14) 



Here we make the explicit computation since it is more difficult than the 
computations for the other first integrals. 

Proof. Let us compute K\ using Th. [3l in this case the parameter A is v 
the first component of the velocity v = (v, 0, 0) which appears in (|2.13j) . In this 
case the computations are easier if we use polar coordinates: we have 



£ = - 
Moreover we have that 



-\Vu\ z + W(u) + [d t S+-\VSf \u 



1 



x„ = x — vt 

t v /, 



and recalling (|2.14|) . we have that the representation T 3v acts on a state ip of 
the system as follows 

(T fl » (*, x)=if>(t,x- vt) e *(v*-*« a *) =u(t,x- vt) e iS(t,x- V t) e i( V -x-iv 2 t) 
namely 

u v (t, x) = u(t,x — vt) 

S v (t, x) = S (t, x — vt) + vxi — 2 v2 t 



Then, by (I2~23D 

PK 1 



dC du v dC dS v ^dt v 
du t dv dS t dv dv 



> ! dS v 
u„ I -q — t + Xl—Vt 
OX\ 



dS v 
dv 



, dS v 

= I Xx - t- | u 

OX\ 



v=0 
2 
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Thus 



Ki ■. ■ / ,r «-'</./• / / — - u 2 dx = j ,r, irtl.r - i P\ 



□ 

The three components of K are the integrals of motion relative to the Galileo 
invariance. Let us interprete this fact in a more meaningful way. If we derive 
both sides of (|4.14p with respect to t, we get 

0=^([ xu 2 dx) - P (4.15) 



dt 

If we define the barycenter (or hylecenter) as follows 

f ~x.u 2 dx f x.u 2 dx 

^j^r-hr- (4 - 16) 

By (f4~Ti5"l) we get 

P 

q = — = const. (4-17) 
7t 

Then, the three integrals of motion relative to the Galileo invariance imply 
that the three components of q are constant. 

From now on, we will study the existence of solitary waves assuming V(x) = 
0. If W S C 2 , we can write 

W(s) = - as 2 + N{s) (4.18) 

where N"(s) = 0. 

Now let us compute the limit (|3.8p . 

Theorem 16. Assume that W is given by j4-18\ ) and that Eq is given by \3.8\) . 
Then 

E = ~a. (4.19) 

Proof. We have that 

em H(\vu\ 2 + \vs\ 2 u 2 )+w^)dx 

a m = y ' = v 

1 ' nm ju 2 dx 

> / \au 2 + N{u)dx _ 1 j N(u)dx 



J u 2 dx 2 J u 2 dx 

Then since N(u) = 0(u 3 ) for u — * 0, we have that 

A m = lim inf > -a. 

In order to prove the opposite inequality, take ^! e .r = sun where ur{x) is 
defined as follows 

1 if \x\ < R 

u R {x) = { if \x\> R+l 

l + R-\x\ if R<\x\ < R+l 
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Then 



da; 



52/ 



u 2 R dx 



J I|V Ufl | 2 + ±. 



dx 



I U R dx 



1 l/|Vu fl | 2 dx / ^iV(ew)dx 
< 2 a+ 2 /u^x +_ R 
1 



1 + 



w R dx 



R 



0{s) 



and so, we have that 



□ 



mi — ; — - < -a 



By gTS} and 1(339]) . we can write 

W(s) = E s 2 + N(s), N(0) = N'{0) = N"(0) 
It is not restrictive to assume that 

E = 

In fact, if ■00 (t, x) is a solution of the equation 



then 

satisfies the equation 



ipi(t,x) = i/j (t,x)e 



-iEnt 



at 2 v 2 



(4.20) 



(4.21) 



(4.22) 



So equation (I4.21|) and (|4.22|) have the same dynamics except of a phase 
factor. This fact, in the swarm model, can be interpreted saying that the internal 
energy of a particle (namely Eq) does not affect the dynamics of the particle. 



4.3 Swarm interpretation of NS 

Before giving the swarm interpretation to equation [NS] we will write it with the 
usual physical constants m and h : 

dtb h 2 1 
^ = -^M+^W'^) + V{x)^ (4.23) 
at 2m 2 

Here m has the dimension of mass and h, the Plank constant, has the dimension 
of action. 
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The polar form of a state tp is written as follows 

i/;(t,x)=u(t,x)e i3 <- t '*V n (4.24) 
and equations (|4. 5|) and (14. 6[) become 

- Au + \w'{u) + (d t S + i | VS| 2 + V(x)] u = (4.25) 
2m 2 \ 2m / 

a t (u 2 )+V- (« 2 ^) =0 (4.26) 

The continuity equation (|3.12j) for INSI is given by (|4.26[) . This equation 

allows us to interprete the matter field to be a fluid composed by particles 
whose density is given by 



PH = u 2 



and which move in the velocity field 



m 

So equation (|4.26p reads 

dtpn + V ■ {pnv) = 

If 



v = ^. (4.27) 



h 2 l 

-— Au+ -W'(u) <^u, (4.28) 
2m 2 

equation (|4. 25(1 can be approximated by 

d t S+ — \\7S\ 2 + V(x) =0. (4.29) 
2m 

This is the Hamilton- Jacobi equation of a particle of mass m in a potential field 
V (cf. eq. I2.44p . The trajectory q(t) of each particle, by (|4.27[) satisfies the 
equation 

q = ™ (4.30) 
m 

(cf. eq. EHD . 

If we do not assume (|4.28p . equation (|4.29[) needs to be replaced by 

d t S+^-\VS\ 2 + V + Q(u) = Q (4.31) 
2m 

with 

-{h 2 /m) Au + W'(u) -(V/m) Au + N'(u) 
Q(u) = = E + (4.32) 

The term Q(u) can be regarded as a field describing a sort of interaction 
between particles. The term Eo does not affect the dynamics, while the other 
terms might contribute to the formations of hylomorphic solitons. 
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Given a wave of the form (|4.24[) . the local frequency and the local wave 
number are defined as follows: 

U)(t, x) = - 



k(t,x) 



H 

VS(t,x) 



n ' 

the energy of each particle moving according to (|4.29|) , is given by 

E = -d t S 



and its momentum is given by 
thus we have that 



P = VS; 

E = Hlu 
p = fik; 

these two equations are the De Broglie relation. It is interesting to observe that 
they have be deduced by the swarm interpretation of the Schroedinger equation. 



We recall again that the swarm interpretation is just a useful pictorial way to 
look at to NS. In physical models, in general, there are different interpretations. 
Here we will mention very shortly some of them. 

In the traditional model of quantum mechanics for one particle, we have 
W = and pn is interpreted as a probability density of the position of this 
particle. 

One of the most important applications of NS is in the Bose-Einstein con- 
densate. In this case W{ip) = ^ H 4 (where U is a constant) and INSI takes 
the name of Gross- Pitaevski equation. Here \ifj\ 2 is interpreted as the particle 
density which in this case are bosons (as, for example, atoms). 



4.4 An existence result for an elliptic equation 

Many existence theorems of solitary waves reduce to the following elliptic equa- 
tion in R N 

-Au + G'(u) = 

u > y J 

where G is a C 1 -function with G(0) = 0. 

This equation has been studied by many authors (see e.g. [57) . [Hi] . 
and their bibliography). In particular, in 14], there are necessary and "almost 
sufficient" conditions for the existence of "finite energy" solutions. 

Since this equation is the basic equation for the existence of solitary waves, 
we give an existence proof. This proof is a variant of the proof in [Mj. It is 
simpler but it uses slightly more restrictive assumptions. 

Equation (|4.33|) is the Euler-Lagrange equation relative to the functional 

J{u) = ^J \X7u\ 2 dx + J G(u)dx, 
We assume that G : M + — > K satisfies the following assumptions: 
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• (G-i) G(0) = G"(0) = 0. 

• (G-ii) G'(s)>c 1 s-c 2 s p - 1 , 2 <p < 2*, s > 0, c x ,c 2 > 0; 

• (G-iii) 3s e M+ : G(s ) < 0. 

Theorem 17. Assume that G satisfies (G-i), (G-ii), (G-iii). Then eg. J^.ggp 
has a nontrivial finite energy solution. 

In order to prove the above theorem, first of all we define an auxiliary func- 
tion G as follows: 
If 

\G'(s)\ <c 3 + c 4 | S r 1 (4.34) 
(where p is defined by (G-ii)), we set 



G(s) 



G (s) for s > 
for s < 



If G does not satisfy (I4.34[) . then, by (G-ii), there exist Si > sq, such that 
G'(si) > 0. In this case we set 

G(s) for 0<s< Sl ; 

G(s) = { G(s 1 ) + G / (s 1 )(s-s 1 ) for s> Sl ; 

for s< 0. 

In any case, G satisfies the assumptions (G-i), (G-ii) (G-iii) and the following 
ones: 



(G-iiii) |G'(s)| <c 3 + c 4 |s|^ 1 , 2 < p < 2* 
(G-v) Vs < 0, G(s) = 



Lemma 18. Let u £ _ff 1 be a solution of the following equation: 

-Au + G'(u) = (4.35) 
Then, u is positive and it is a is a solution of l{4-33\ ). 

Proof. The fact that u is positive is a straightforward consequence of the 
maximum principle and the fact that G'(u) = for u < 0. Then if (|4.34D holds, 
G'(u) = G'(u) and hence u is a solution of (|4.33|) . If (|4.34p does not hold, 
by the maximum principle it follows that u{x) < sx; then, also in this case 
G'(it) = G'(u) and u is a solution of (|4.33p . 

□ 

Thus, by the above lemma it is not restrictive to assume that G satisfies 
(G-i), (G-iii), (G-iiii) and (G-v), since otherwise we can work with G. 
We set 

H^lueH 1 (R N ) : u = u(|a?|)> 
Lemma 19. There exists u £ fL\ such that 

J G{u)dx < 
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Proof. Wc set 



so for \x\ < R; 

ur(x) = { s -s (\x\-R) for R<\x\<R+l; 
for \x\>R+l. 



Thus we have 



G(u)dx 



< 



R+l f R 



<C X 



max G(s) / r N - x dr+ I G(s )r"- L dr 
\s\G[R,R+l] W J J R Jo 

(R+lf-R N ]+±G(s )R N 



1 



(4.36) 



< C 2 R N - L + -G(s )R 

where C\,Ci are positive constant. By (G-iii), .G(sn) < 0; hence, for R suffi- 
ciently large, J G(un)dx < 
□ 

Lemma 20. If u e 2? EL 

is a radially symmetric function, then 
\u(x)\ < C N — — 



Proof. We have that 
d 



(r^u 2 ) = 2r N - x u^ + (N - 1) r"" V > 2r Ar " 1 u^ 
fir v ' dr ar 



Then, integrating over (R, +oo) we get 

-R N - 1 u(R) 2 > 2 / 



r u—dr 
dr 



and so 



/•+0O 

R N - 1 u{R) 2 < 2 / 



<ir 



- N ~ l dr < 



-t-oo 



du 
dr 



+ u 2 )r N - 1 dr<C N \\u\\ z Hl 



□ 



Theorem 21. (Strauss For p S (2,2*) i/ie embedding 

Hi L p 

is compact. 



3 Here I? denotes the space of C°°-functions with compact suport. 
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Proof. Let u n — 1 weakly in if*; we need to prove that u n — ► strong 
in LP . Since u n 0, then there is a constant M such that ||u„|| ff i < M; so by 
lemma [501 we have that 



, IIP" 2 / Li I 2 



_j / \u n \ 2 dx (4.37) 

( J ^fe-V / Kl^ (4.38) 

M V~ 2 „ 1|2 . AfP 



where 



(JV-l)(p-2) >Q 



2 

Since u„ — ► strongly in L p (Br) , by (|4.37|) . we have that 

liminf / \u n \ p dx= limmf\\u n \\ p LP(B ■. + liminf||u„||^ „_ , < 

n— too J n— too v^ni n—>oo ^ "RJ 

By the arbitrariness of R, it follows that this limit is 0. 
□ 

Proof of Th. [T71 Take a function f3 G C X (R) such that 

• -(/3-i) Vs G K, < /3(s) < 1 

• Vs < 0, 0(s) > 

• -03-iii) Vs > 0, /3(s) = 

• -(/3-iiii) Vs G M, /3'(s) < 0. . 

Now we define a C 1 -functional on Hi as follows: 



MP 



F{u) = \ J |Vu| 2 dx - bp (J 



G(u)dx I ; 



here & is a positive constant defined by 



y\\/u\ 2 dx + l 



13 (J G{u)dx) 

and u is defined by lemma [TO1 This choice of b, implies that 

F(u) = -1 

Clearly F is a C 1 functional defined on if*. 

Since /? is bounded, then F(u) is bounded below and its infimum is a num- 
ber less or equal to —1. Then there exists a minimizing sequence u n . Also, 
bounded implies \ J |Vu n | dx bounded and hence u n — > w weakly in T>, 



so, 3 



A T>r' 2 denotes the closure of the set of radially simmetric C°°-functions with compact 
support with respect to the norm </ J | ^Ju\ 2 dx. 
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Since the infimum is negative we have that f3 (J G{u n )dx) < 0, and hence, by 



G(u n )dx < 0. (4.40) 
By (G-ii), we have that there is a constant C5 > such that, for s > 

G(s)> f G'(t)dt> f (at - c 2 t p - 1 ) dt 



jo 



^1 2 c 2 „ l 2 2* 

> -CiS S F > -C4S — C5S 

- 2 p ~ 2 

By this inequality and (|4.40|) . we get 

0> y G(u n )dx > -c 4 J\u n \ 2 dx-c 5 J \u n \ 2 dx 

> \ c ± IWIz,= - c e llVunll^j 

Thus since || Vu„|| i2 is bounded, by the above inequality, also ||u n || i 2 is bounded 
and hence also ||it n || ffl is bounded. So, by theorem 1211 u n — > u> in L p strong, 
where p is defined by (G-ii). By our assumptions, the functional 



G(u)dx 

is continuous in L p . Then we have that J G(u n )dx — ► J G(w)dx and hence 
/3 (J" G(u n )dx) — ► /3 (/ G{w)dx) . Since ^ J |Vu| 2 (ix is l.s.c. it follows that u> 
is the minimum of _F. Thus w is different from and it satisfies the following 
equation: 

-Aw + XG'(w) = 

where 



G(w)dx 



By (/3-iiii), we have that A > 0, and since w 7^ 0, we have that A > 0. Now set 

«(*)=«(^) 

Clearly u satisfies equation (|4.33|) 
□ 

4.5 Standing waves and travelling waves 

In any complex valued field theory, the simplest possible solitary waves are the 
stationary waves, namely finite energy solution having the following form 

4(t,i) = u(i)r iu ',u >0, (4.41) 

In particular, for the nonlinear Scroedinger equation (|NS[) . substituting 
in eq. (|NS[) . we get 

- Au + W'(u) = 2uu (4.42) 
Now we can apply the results of section 14.41 and get the following theorem: 



35 



Theorem 22. Assume that W has the form $4.20\ ) where N satisfies the fol- 
lowing assumptions: 

N'(s) > -cs^ 1 for any s > 1, (4.43) 

3s G R + : N(s ) < 0. (4.44) 

then eq. tNb}) has finite energy solitary waves of the form ft4-41\ ) for every 
frequency lo G (Ei,Eq) where Eq is defined by {4-20^ and 

Ei = inf {a G R : 3s G R + , as 2 > E s 2 + N(s)} 

Notice that, by virtue of (|4.44[l . E\ < Eq. Also, it is possible that E\ — —oo; 
for example this happens if 

Wis) = — \s\ p , 2<p<2* (4.45) 
P 

Proof. We want to apply theorem [T71 to eq. (|4.42[k to this end we set 

G(s) = W{s) - us 2 = {E Q -lj)s 2 + N(s) 
It is trivial to verify that, for every tu G (Ei,Eq) , G satisfies the assumptions 

of Th. mi 

□ 

Now we will exploit the other symmetries of equation INSI to produce other 
solutions. First of all, since INSl is invariant for translations, for any xq G M. n , 
the function 

ip X0 {t,x) =u(x-x )e-^ t (4.46) 

is a standing wave concentrated around the point xq. The space rotations do 
not produce other solutions since u(x) is rotationally invariant. 

Since the Lagrangian (|4. 3[) is invariant for the Galileo group, we can obtain 
other solutions: we can produce travelling waves just applying the transforma- 
tion (pTHf to (j4~l6l) 

lk ,v(t,a:) = u{x - x - vty (v - x ~ Et \ E = ^v 2 +lu 

ipxcvitjx) is a solitary wave concentrated in the point xq + vi, an hence it 
travels with velocity v. 

Finally, other solutions can be produced by the invariance (|2.15| : for 9 G 
[0, 2tt) , we have the solutions 

lko,v,fl(f, x) = u(x -xq- vty^ x - Et+ ^ (4.47) 

The invariance by time translations do not produce new solutions, since a time 
translation on ijjx ,v.e produces a space and phase translation. Concluding, if 
we fix a charge a, we obtain a radially symmetric solution of the form (|4.41jl : 
by the invariance of the equation, this solution produces a 7-parameters family 
of solutions given by 14.471 

If we consider lNSI with the usual physical constants m and h, namely equation 
(|4~2"3"1) with V = 0, (TP?]) takes the form 

tfx 0l v,fl(t, x) = u(x -xq- vt) e *(p-*-B*)/«+« (448) 



3G 



where 



p = mv 

1 , 
E = —mv + ujH 

The meaning of these relations within the swarm interpretation is the following: 
^x ,v,e{t, x) is interpreted as a swarm of particles of mass m; p is the momentum 
of each particle and E is its energy: ^mv 2 is its kinetic energy and u>K is related 
to the nonlinear interaction of the particles. 



4.6 Existence of solitons 

Theorem 1221 provides sufficient conditions for the existence of solitary waves. In 
order to prove the existence of hylomorphic solitons (cf. Def . [TT)|) , first of all it 
is necessary to assume that the Cauchy problem for INSI is well posed namely, 
that it has a unique global solution which depends continuously on the initial 
data. For example, this is the case if W is a globally Liptschitz function; we 
refer to the books [55] and [T7] for more general conditions. 

Moreover, it is necessary to investigate under which assumptions the energy 

£ = J (~IW>| 2 + wW) dx 

= J Q|Vu|Vi|VS|V + T4^))d* 



achieves the minimum on the manifold 



\ip\ 2 dx = a 



Clearly, if u minimizes the functional 

1 



J(u) = j ^\Vu\ 2 + N{u)dx 



M rT = {u£ H l (R N ) : / \u\ 2 dx = a 



on the manifold 



then the set of all minimizers of 8 is given by ip{ x ) = u(x)e lB , 9 € [0, 2ir) 
We make the following assumptions on iV : 

\N'{s)\ < cxls^- 1 + c 2 \s\ p - 1 for some 2 < q < p < 2*. (4.49) 

N(s) > -cis 2 - c 2 |s| 7 for some c u c 2 > 0, 7 < 2 + — (4.50) 

Theorem 23. Let N satisfy U4^ , U^Uj) and [f44\ ). Then, 3 a such that, 
V <7 > a, J has a minimizer u on M a which is positive and radially symmetric 
around some point. Moreover, if the Cauchy problem is well posed ip := ue is 
a hylomorphic soliton. 
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In order to have stronger results, we can replace (|4.44[) with the following 
hypothesis 

N(s) < -s , 2 < [3 < 2 + — for small s. (4.51) 

In this case we find the following results concerning the existence of the 
minimizer of J(u) for any a. 

Theorem 24. // |^.^9[ ) ; {4-50fy and |^.5J[ ) hold, then the same conclusions of 
Th\M hold for every a > 0. 

In particular, for N — 3 we have 

Corollary 25. Let N = 3. If (44<fy and (4~50\ ) hold and N g C 3 , with N'" (0) < 
0, then the same conclusions of Th \23\ hold for every a > 0. 



For the proofs of ThUjH Thailand ThUHl we refer to jlj. Analogous results 
have been obtained by Cazenave and Lions [TH] when W(u) = ^\u\ p . 

4.7 Dynamics of solitons 

In this section we will describe a recent result relative to the dynamics of soli- 
tons when a potential V(x) is present. As we will see, we will introduce three 
parameter, h, a and 7 in our equation and we obtain a meaningful result for 
particular values of these parameters. First of all, let us consider the following 
"unperturbed" Cauchy problem: 

lh lti = + 2^ w ' {h ^ ] (4 ' 52) 

^(0 iX ) = ±u(^yt™ (4.53) 

where 

/^l + ^y 1 (4-54) 

and U : M. N — > M, N > 2, is a positive, radially symmetric solution of the 
equation 

- AU + W'(U) =2uU (4.55) 

with 

We set 



\\U\\ L2 =a (4.56) 



u h (x) = h^U (p) 



and establish a relation between a, (3 and 7 in order to have stationary solution 
of (|4.52|) of the form ip(t,x) = U} l {x)e^ t ~^ t : replacing this expression in (|4.52j) 
we get 

- h 2 Au h + ±-W'(h~<u h ) = 2cj h u h . (4.57) 
If we take the explicit expression of u%, we get 
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and hence, by rescaling the variable x, 

-tf-i-W+<*MJ{x) + W (U(x)) = 2oj h h a ^U (x) . 
Thus, comparing the above expression with (|4.55p . we get (|4.54[) and 

uj h = (4.58) 

Using the arguments of section [431 we have that the solution of (|4.52[) . (|4.53|) 
is given by 

titrf-^uf '-*-* ) ei(— «) (4.59) 



hi V hP 
with 

T? 1 2 ■ W 

h = -v 



2 h a -~< 

Moreover if the problem (j4.52|) . (|4.53|) is well posed this is the unique solution. 

We can interpret this result saying that the barycenter q(t) of the solution 
of (g3H), (14*35)) (defined by [OH) satisfies the Cauchy problem 

q = 
9(0) = 9o 
9(0) = v 

Let us see what happens if the problem (|4.52p , (14.53[) is perturbed, namely 
let us investigate the problem 



dib h 2 1 



if} (0,x) = ifh(x) 



where 

and wo is small, namely there is a constant C such that 

WMrn < Ch a ~^, 
V(x)\w (x)\ 2 dx < Ch a -i. 



(4.60) 



(4.61) 



Also we assume that 

\\U + w \\ L2 = \\U\\ L2 =a 

Theorem 26. Let all the assumptions of Th. \23\ hold and let a > a (a is 
defined by J^.56'[ ) ). Moreover we assume that V : V : R N -> K is a C 2 function 
such that: 

V(x) > 0; (V ) 
\W(x)\ < V(x) b for \x\ >R 1 >l,b£ (0, 1); {V x ) 
V(x) > \x\ a for \x\ >R 1 >l,a>l. (7 2 ) 
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Finally, we assume that 

a > 7 (4.62) 

Then the barycenter (t) of the solution of the problem J^.6'fl[ ) satisfies the 
following Cauchy problem: 

q h (t)+VV(q h (t))=H h (t) 
Qh (o) = q 
<7fc(0) = v 

where 

sup\H h (t)\ -> as 

Proof: The proof of this theorem can be found in [12] . 
□ 



Remark 27. The assumption ( IV2D is necessary if we want to identify the po- 
sition of the soliton with the barycenter ()4.16|) . Let us see why. Consider a 
soliton ip(x) and a perturbation 

ij, d {x) = ip(x) +(p(x-d), d e R N 

Even if ip (x) -C ip{x), when d is very large, the "position" of ip(x) and the 
barycenter of ipd(x) are far from each other. In [T^] (lemma 25), it has been 
proved that this situation cannot occur provided that QV2D hold. 

Remark 28. We will give a rough explanation of the meaning of the assumption 
a > 7 which, in this approach to the problem, is crucial. The energy, using eq. 
(|4.8[) can be divided in two components: the internal energy 

J h {u)= J (y|Vu| 2 + W h (u))dx (4.63) 

and the dynamical energy 

G(u,S)= J (^\VS\ 2 + V{x)^ju 2 dx (4.64) 

which is composed by the kinetic energy | J |V5| 2 u 2 dx and the potential energy 
J V(x)u 2 dx. By our assumptions, the internal energy is bounded from below 
and the dynamical energy is positive. As h — > 0, we have that 

Jh (V/0 = h N P- a -^ 

and 

Then, we have that 

Jh (iph) 

So the assumption a — 7 > implies that, for h <C 1, G (iph) <C (V'fc); namely 
the internal energy is bigger than the dynamical energy. This is the fact that 
guarantees the existence and the stability of the travelling soliton for any time. 
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We end this section with an heuristic proof of ThJ^U This proof is not at 
all rigorous, but it helps to understand the underlying Dynamics. As in section 
we interpret p-u = u 2 as the density of particles; then 

TL = I pndx 



is the total number of particles. By (|4.3ip . each of these particle moves as a 
classical particle of mass m — 1 and hence, we can apply to the laws of classical 
dynamics. In particular the center of mass defined in (|4.16|) takes the following 
form: 

= Jxrnpndx = Jxpndx_ 
J mpudx J pndx 

The motion of the barycenter is not affected by the interaction between particles 
(namely by the term (I4.32[> ) . but only by the external forces, namely by VV. 
The global external force acting on the swarm of particles is given by 



F = - J VV{x)p H dx. (4.66) 

Thus the motion of the center of mass q follows the Newton law 

~F = Mq, (4.67) 

where M = J mpudx is the total mass of the swarm; thus by (|4. 65[) . (|4.66|) and 
PUT)) , we get 

_ JVVpndx _ fVVu 2 dx 
m J pndx m J u 2 dx 

If we assume that the u(t, x) and hence pn,(t, x) is concentrated in the point 
q(t), we have that 

J VVu 2 dx = VV (q(t)) J u 2 dx 

and so, we get 

mq(t) = -VV (q(t)) . 

Notice that the equation mq(t) = — W (q(t)) is the Newtonian form of the 
Hamilton- Jacobi equation (|4.29[) . 



5 The nonlinear Klein-Gordon equation 

5.1 General features of NKG 

The D'Alambert equation, 

□V> = o 

is the simplest equation invariant for the Poincare group, moreover it is invariant 
for the "gauge" tranformation 

Tp 1— ► ijj + C 

Also, ifip is complex valued, it is invariant for the action (|2.15p . Thus, it satisfy 
assumptions A-l, A-2 and A-3, but it is linear and it does not produce solitary 
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waves. There exist only non-dispersive waves. Let us add to (|2.12[) a nonlinear 
term: 

C= l -\d t M 2 ~\WM 2 -W(^>) (5.1) 

where 

W : C -> R 

satisfies the following assumption, 

W" (e^V) = ^ {tp) 

namely W (tp) = F(\ip\) for some function F = K — > K. This is simplest non- 
linear Lagrangian invariant for the Poincare group and the trivial gauge action 

The equation of motion relative to the Lagrangian (|5.1[) is the following: 



Dip + W (ip) = (NKG) 



where W (tp) is intended as in 

In the following sections we will see that equation INKG) with suitable (but 
very general) assumptions on W produces a very rich model in which there are 
solitary waves and solitons. Moreover we will see that these solitons behave as 
relativistic particles. 

If W'(ip) is linear, namely W'(ip) — m 2 ip, then eq. INKGI reduces to the 
Klein-Gordon equation 

aip + m 2 ip = 0. (5.2) 

Among the solutions of the Klein-Gordon equations there are the wave packets 
which behave as solitary waves but disperse in space as time goes on. On the 
contrary, if W has a nonlinear suitable component, the wave packets do not 
disperse and give hylomorphic solitons. 

Sometimes, it will be useful to write ip in polar form, namely 

ip{t,x) =u(t,x)e iS ^ x) (5.3) 
In this case the action J Cdxdt takes the fom 

S(u,S) = ^J (d t u 2 ) - \Vu\ 2 + (d t S 2 ) - |VSf u 2 dxdt - J W(u)dxdt = 

(5.4) 

and equation INKGI becomes: 



Uu- 



(d t S 2 ) - \VS\ 2 u 2 + W'(u) = (5.5) 
d t (u 2 d t S) - V • (u 2 WS) = (5.6) 
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5.2 First integrals of NKG and the hylenic ratio 

Easy computations and the results of section 12. 71 show that the integral of 
motion of lNKGI are given by the following expression: 



• Energy. We get 

t 

Using (|5.3p we get: 
£ -- 



\ (dtuf + \ |V,i| 2 + \ [(d t S) 2 + \VS\ 2 ] u 2 + W(u) 



(5.7) 



dx (5.8) 



• Momentum. We have 



P = - Re / <9 t V>V-0 dx 



Using (|5.3[) wc get: 



J (dtuVu + dtSVS u 2 ) dx 



• Angular momentum. We have 

M = Re / xx Vipd^j) dx 



Using (|5.3p we get: 

M = J (x x VS* d t S u 2 + x x Vu d t u) dx 

• Hylenic Charge. We have 

7i = Im / dtipip dx 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



Using (|5.3p we get: 



H = / d t Su 2 dx 



(5.13) 



• Ergocenter velocity. If we take the Lagrangian ()5.1|) . the quantity 
preserved by the Lorentz transformation, is the following 



K = tP - / x 



\\dM 2 + \ |W| 2 + W(il>) 



dx (5.14) 



43 



The computation of K is more involved than the prrevious ones and we will 
make it in details: 

Proof. Let us compute Ki using Th. [3J in this case the parameter A is the 
velocity Vi which appears in ()2.8|) : we have 



dC gV 
dipt dvi 



p Ki = Re 

= Re I dti> 



C 



dv. 



dtp dt ^ ® x k 
~dldv~ + ^ 



k=l 



dx k dv l 



where the derivative with respect to Vi need to be computed for Vi = 0. 



Since for k ^ i, = 0, we have that 



dt_ 

dvi 



Re (d t ipd Xi ip) 



dxj 
dv. 



\\dti>\ 2 + \\vi>\ 2 + w(*p) 



= pe-R 

ovi 

Also we have 



PPi ■ 



dv. 



Re (dti>d Xt ip) 



dv 



dvi 



dt_ 

dvi 



dt_ 

dvi 



1 / Vi=0 



dx 
dvi 



1 / Vi=0 



so that 



d t 



dvi ^1 - vf 
d x — vd 



PK t 



dvi y/\ - v f 



-P£Xi + PP,t 



Integrating in the space we get the «-th component of (|5. 14[) . 

□ 



Let us interpret (j5. 14|) in a more meaningful way. If we derive the terms of 
()5.14j) with respect to t, we get 



dt 



l\d^\ 2 + \\V^\ 2 + WW 



dx = I xpgdx 



Now, we define the ergocenter (or barycenter) as follows 

Jx [i \dM 2 + \ |v^| 2 + ww] dx j xpsdx 



Q 



dj- 



then, by the conservation of £ and eq. (|5.15[) . we get 

«4 



(5.15) 



(5.16) 



(5.17) 
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Concluding, the Poincarc group provides 10 independent integral of motions 
which are £, P, M, K; they can be replaced by integral of motions £, P, M, Q 
since also these quantities are independent. 

Notice the difference between (|4.14p and (|5.14p and consequently the dif- 
ference between the hylecenter (|4.16j) and the ergocenter (|5.17|) . In INSI the 
barycenter or "center of mass" or coincide with the hylecenter or "center of 
hylenic charge" ; in INKG1 the barycenter coincide with the ergocenter or " cen- 
ter of energy" . 

The precise definition of mass and its meaning will be discussed at pag. 1511 



We now assume that W is of class C 2 and we set 
W(s) = Kn 2 s 2 + N{s) 

where m 2 = W"(0). 

Theorem 29. IfW is of class C 2 , then we have that 



(5.18) 



E 



lim inf 



s^o *ex £ |W(*)| 

Proof. We have W = (tp^t) = (u,u t> w,k) ; then by JSH), ([5~T3]) and (??) 



> 



I 


\ (d t uf 


+ \ \Vu\ 2 + \ [(d t S) 2 + (VS) 2 ] u 2 + W(u) 


dx 


I [\^n 2 + 


\J ui u 2 dx\ 
\m 2 u 2 + N(u)] dx 





J \u>\ u 2 dx 



Since 



|o;| u dx < uj u dx 



1/2 



L 2 dx 



1/2 



< 



1 

m 

1 
2m 

1 

2m 



j 2 u 2 dx 



1/2 



i 2 u 2 dx 



1/2 



ui 2 u 2 dx + / m 2 u 2 dx 



- — / (ui 2 + m 2 ) u 2 dx 
2m J ' 



we have that 



S (*) / [\uj 2 u 2 + \m 2 u 2 + N(u)] dx 



J N(u)dx 



\H{9)\ '- ^j(u; 2 + m 2 )u 2 dx 
Then since N(u) = 0(u 3 ) for u — > 0, we have that 

lim inf - — - — \- > m 

£ ^0 4. e x e |H(*)| ~ 



tt— f (oj 2 + m 2 ) u 2 dx 
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In order to prove the opposite inequality, take ^s.r = (Sun,—iSuR) 
(u, 0, 1, 0) where 



1 if \x\ < R 

ur(x) = { if \x\ > R+ 1 

1 + R- \x\ if R < \x\ < R+ 1 



(5.19) 



Then 



g(g) £ (j^g) _ £2 I[W u r\ +\ur + ^W(eu) 



da; 



d:r 



= 1 + 



R 



f u 2 R dx 
f 0(e) 



< 



l J\X7u R \ 2 dx f^N(eu)d 



2 J u 2 R dx 



I u R dx 



□ 



5.3 Swarm interpretation of NKG 

Before giving the swarm interpretation to equation INKG1 we will write it with 
the usual physical constants c, m and h : 



with 



d 2 j> 



W{u) 



c 2 &ip + W'(ip) = 



■u 2 + N(u) 



(5.20) 



2h 2 



Here c has the dimension of a velocity (and represents the speed of light), m 
has the dimension of mass and H is the Plank constant. 
The polar form of tp is written as follows 



i(j(t,x) = u(t 7 x)e zS ^ x)/R 
and equations (|5.5p and (|5.6p become 



(5.21) 



(d 2 u - c 2 Au + N'(u)) + (-<9 t 2 S* + c 2 |V5| 2 + m 2 c 4 ) u = (5.22) 



d t (u 2 d t S) - c 2 V • (u 2 VS*) = (5.23) 

The continuity equation (|3.12[) for INKGI is given by (|5.23|) . This equation 
allows us to interprete the matter field to be a fluid composed by particles whose 
density is given by 

p H = ~u 2 d t S 



and which move in the velocity field 



tPdtS' 



(5.24) 
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If 

h 2 (d 2 u - c 2 Au + N'(u)) < (-<9 t 2 S + c 2 |VS| 2 + m 2 c 4 ) u, (5.25) 

namely, if fi is very small with respect to the other quantities involved, equation 
(|5.22p can be approximated by 

d 2 S = c 2 \VS\ 2 + m 2 c i . (5.26) 

or 



d t S+\Jm 2 c A + c 2 \VS\ 2 = (5.27) 

This is the Hamilton-Jacobi equation of a free relativistic particle of rest mass 
m (cf. eq. I2.50[) whose trajectory q(t), by (|5.24[) satisfies the equation 

V <? 

« - "355 < 5 ' 28 > 

(cf. eq. [23B 

If we do not assume (|4.28|) . equation ()4.29[> needs to be replaced by 



dtS = ±\jm 2 c A + c 2 \VS\ 2 + Q(u) 

with 

Q(u) = fi 2 (<9 t 2 u - c 2 Au + 7V'(u)) 

The term Q(u) can be regarded as a field describing a sort of interaction 
between particles. 

Given a wave of the form (|5.21[) . the local frequency and the local wave 
number are defined as follows: 

d t S(t,x) 
u>(t,x) = — - 



k(i, x) 



h 

VS(t,x) 



the energy of each particle moving according to (|5.27p , is given by 

E = d t S 



and its momentum is given by 
thus we have that 



P = VS; 

E = hw 
p = fik; 

these two equations are the De Broglie relation. 
Thus, eq. (|5.24[) becomes 



v = - = f (5.29) 
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5.4 Existence of solitary waves and solitons in NKG 

The easiest way to produce solitary waves of INKGI consists in solving the static 
equation 

-Aii + W'(u)=0 (5.30) 



and setting 



(X\ — vt 
2 ,X2,x 3 ) ; (5.31) 



ip v (t, x) is a solution of cq. INKGI which represents a bump which travels in the 
cci-direction with speed v. 

Thus by Theorem \T7\ we obtain the following result: 

Theorem 30. Assume that W satisfies (G-i), (G-ii), (G-iii). Then eq. \NKG\ 
has real valued solitary waves of the form \5.31}) . 

However, it would be interesting to assume 

W > 0; (5.32) 
in fact the energy of a solution of equation INKGI is given by 



l\d^\ 2 + ±\\7^\ 2 + W^) 



dx 



Thus, (|5.32|) implies that every state ip has positive energy. In this case, the 
positivity of the energy, not only is an important request for the physical models 
related to this equation, but it provides good a priori estimates for the solutions 
of the relative Cauchy problem. These estimate allows to prove the existence and 
the well-posedness results under very general assumptions on W. Unfortunately 
Derrick [35] , in a very well known paper, has proved that request (|5.32|) implies 
that equation (|5.30|) has only the trivial solution. His proof is based on the 
following equality (which in a different form was also found by Pohozaev ). 
The Derrick-Pohozaev identity (see e.g. ([7])) states that for any finite energy 
solution u of eq. INKGI it holds 



1 1 

N ~ 2 



J \Vu\ 2 dx + J W{u)dx = (5.33) 



Clearly the above inequality and (|5.32|) imply that u = for N > 2. 

However, we can try to prove the existence of solitons of eq. INKGI (with as- 
sumption (|5.32p ) exploiting the possible existence of standing waves (as defined 
by (|4.41[) ), since this fact is not prevented by ea. (|5.33[) . 
Substituting (|4~4"Tj) in eq. INKGI we get 

- Au + W'(u) =lu 2 u (5.34) 

Since the Lagrangian (|5.ip is invariant for the Lorentz group, we can obtain 
other solutions ipi(t,x) just making a Lorentz transformation on it. Namely, if 
we take the velocity v = (v, 0, 0), |u| < 1, and set 

t' = 7 (t — vxi) , x\ = 7 (x\ — vt) , x' 2 = X2, x' 3 = X3 with 7 



VT 
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it turns out that i[>i(t, x) = ip(t', x') is a solution of INKG1 

In particular given a standing wave ip(t, x) = u{x)e~ lu)t , the function ip v (t, x) := 
ip(t',x') is a solitary wave which travels with velocity v. Thus, if u(x) = 
u(xi,X2, Xs) is any solution of Eq. (|5.34|) . then 

ipv(t,xi,x 2 ,x 3 ) = u{j(x 1 -vt),x 2 ,x 3 )e l{k - x - UJt \ (5.35) 
is a solution of Eq. INKGI provided that 

u> = 7W0 and k = jujqv (5.36) 
Notice that (|5.3ip is a particular case of (|5.35[) when ujq = 0- 

We write W in the form (|5.18[) and we make the following assumptions: 

• (W-i) (Positivity) W(s) > 

• (W-ii) (Nondegeneracy) W = W(s) ( s > 0) is C 2 near the origin with 
W(0) = W'(0) = 0; W"(0) = to 2 > 

• (W-iii) (Hylomorphy) 3s : N(s ) < 

• (W-iiii) (Growth) there is a constant c > such that 

N'(s)>-c 1 s-c 2 s p - 1 ) 2<p<2* 

Here there are some comments on assumptions (W-i), (W-ii), (W-iii). 

(W-i) implies that the energy is positive; if this condition does not hold, it 
is possible to have solitary waves, but not hylomorphic waves (cf. Proposition 
16 of [ID]). 

(W-ii) In order to have solitary waves it is necessary to have W"(0) > 0. 
There are some results also when W"(0) = (null-mass case, see e.g. [2] and 
[2]), however the most interesting situation occurs when W"(0) > 0. 

(W-iii) This is the crucial assumption which characterizes the potentials 
which might produce hylomorphic solitons. This assumption permits to have 
states $ with hylomorphy ratio A (JS) < m. Actually, N is the nonlinear term 
which, when it is negative, produces a attractive "force". 

We have the following result: 

Theorem 31. Assume that (W-i),..., (W-iiii) hold, then eq. \NKG\ has finite 
energy solitary waves of the form ip(t,x) — u{x)e~' luJt for every frequency lu G 
(too, to) where 

too = inf | a G K : 3u G K + , ^aV > W(u)\ 

Notice that by (W-iii), too < to, then the interval (too, to) is not empty. 

Proof. By the previous discussion, it is sufficient to show that equation 
(|5.34|) has a solution u with finite energy. The solutions of finite energy of 
(|5.34|) are the critical points in the Sobolev space H 1 (R 3 ) of the reduced action 
functional: 

J(u) = i 1 1 Vu| 2 dx + { G{u)dx, G (u) = W{u) - iwV (5.37) 
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Now we apply theorem 1171 

It is easy to check that for every frequency to € (mo, to) , the required as- 
sumptions are satisfied. 

□ 

In [1] the existence of soliton is proved. The proof is quite involved and will 
not discussed here; we only refer that it is necessary to strengthen assumption 
(W-iiii) with the following one: 

• (W-iiii') (Growth condition) Al least one of the following assumptions 
holds: 

— (a) there are constants a, b > 0, 2 < p < 2N / (N — 2) such that for 
any s > : 

\N'(s)\ ^asP- 1 +bs 2 -l. 

- (b) 3si > s a : N'(si) > 0. 

5.5 Dynamical properties 

In this section we will show that the solitons and the solitary waves relative to 
cq. INKGI behave as relativistic bodies. In fact, the relativistic effects like the 
space contraction, the time dilation and the equality between mass and energy 
are consequences of of the variational principle A-l and the invariance for the 
Poincare group. 

First of all observe that, by Eq. (|5.35[) . the following theorem follows: 

Theorem 32. Any moving solitary wave experiences a contraction in the di- 
rection of its movement of a factor 1 /*y with 7 = p . 

Thus the space contraction is a trivial fact. On the contrary, the time dilation 
needs a more subtle computation. 

The standing waves of eq. can be considered as a clock. Let us denote by 
q (t) the position of our clock at the time t. 

If we assume that at t = Q, q(0) = (0, 0, 0), the motion of the clock is given 

by 

q(t) = K0,0); 

then the behavior of the moving clock at the point q(<) is obtained replacing x 
by q(i) in Eq. ^M- 

^(i,q(t)) = V„(Mt,0,0) 

= u (0,0,0) e i( k -q(*)-^*) ; 

taking into account eq. (|5.36[) . we get 

k • q(£) — cot = 'JluqVi ■ V\t — ut 
= ("fujovl - 7^ )i 
= 7^o (vj - 1) t 

= ^t 

7 
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Then 

^„(t,q(t)) =u (0,0,0) e"*^' 
this equation shows that our moving clock is vibrating with a frequency 

W ° /l 2 

7 

Since the intervals of time measured by a clock are inversely proportional to the 
frequency of the vibrations, we get that 

Then we get the following 

Theorem 33. A moving clock moves slower than a resting clock by a factor 
I- 1 - 

In classical mechanics the mass m can be defined as a quantity which relates 
the momentum P = (Pi, P2, P3) to the velocity v = (v±, «2, ^3) by the following 
formula 

P = rnv 

Since the momentum of a solitary wave is defined by (|2.32[) . it is possible to 
define the mass of a solitary wave by the above formula and to compute it. 

In the case of the Schroedinger equation, the velocty of a soliton is given by 
q where q is defined by (|4.16|k then by (|4.17p we have that 

P 

mass = — = H (5.38) 

namely the mass of the soliton equals the hylenic charge. Actually, if we add 
the usual constant as in eq. (14.23p . we get that 

mass = mTL 

This equation is consistent with the swarm intepretation of section 14.31 it says 
that the soliton consists of a number Ti. of particles of mass m. 

In the case of the equation lNKG) the velocty of a soliton is given by Q where 
Q is defined by ()5.16|) : then by (15.17|) we have that 

P , 

mass = — = c 

Q 

namely the mass of the soliton equals its energy. If we add the usual constant 
as in eq. (|5.20|) . we get that 

£ 

mass = — 
c z 

namely the celebrated Einstein equation. Namely, we have the following: 

Theorem 34. The mass of a solitary wave of eq. \NKG\ is proportional to its 
energy with the factor of proportionality cr 2 . 
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